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We study the Gutzwiller method for the spinless fermion model in one dimension, which is one of 
the simplest models that incorporates the intersite Coulomb interaction. The Gutzwiller solution of 
this model has been studied in the literature but with differing results. We obtain the Gutzwiller 
solution of the problem by a careful enumeration of the many-particle configurations and explain 
the origin of the discrepancy in the literature to be due to the neglect of the correlation between 
the neighboring bond occupancies. The correct implementation of the Gutzwiller approach yields 
results different from the slave-boson mean-field theory, unlike for the Hubbard model with on-site 
interaction, where both methods are known to be equivalent. The slave-boson and the Gutzwiller 
results are compared to the exact solution, available for the half-filled case, and to the numerical 
exact diagonalization results for the case of general filling. 
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Introduction - Conventional many-body theory has dif- 
ficulty in taking into account the local correlations while, 
at the same time, preserving the itinerant character of 
the many-electron states. Because of this, there has 
been a resurgence of interest in the variational Gutzwiller 
methodfi where local correlations show up as a renor- 
malization of the kinetic energy^r^ within the Gutzwiller 
approximation and the overall itinerant character of the 
electrons is preserved. The method leads to a many- 
body modification of the Fermi surfaces and is being 
increasingly used to describe the measured Fermi sur- 
face topology of correlated systems starting from ma- 
terials as simple as Ni^ to more complex systems such 
as NaCo02® and LaOFeAs^ using model Hamiltonians. 
This in turn has led to a flurry of work for incorporating 
the Gutzwiller approach in the density-functional theory 
for realistic solids 4"— 

While the Gutzwiller method has been well developed 
for on-site Coulomb interactions, it has been less stud- 
ied when the inter-site Coulomb interactions are present. 
On the other hand, the nearest-neighbor Coulomb inter- 
action is an important ingredient in models used to de- 
scribe the charge order-disorder transition in many sys- 
tems such as the Verwey transition in the half-metallic 
Fe304. Thus there is a need to develop the Gutzwiller 
approach for incorporating the effects of the inter-site 
Coulomb interactions, which presents nuances not en- 
counterd in treating models with on-site Coulomb inter- 
actions such as the Hubbard model. 

The simplest model that contains the intersite 
Coulomb interaction is the spinless fermion model in one 
dimension (ID) 



H 



h.c. 



where c\ creates a spinless fermion at the site i, hi 



(1) 



and the summation is taken over distinct pairs of near- 
est neighbors. The is also an important model due 
to its equivalence to the XXZ Heisenberg model via 
the Jordan- Wigner transformation The Gutzwiller ap- 
proach for this model has been treated by three different 
papers in the literaturei^^— including two by the same 
author; however the results differ from one another in- 
cluding the all-important Gutzwiller reduction factor for 
the kinetic energy. In this paper, we study this model 
with the Gutzwiller approach by a careful enumeration 
of the many-particle configurations. Our results agree 
with one^'^ of the two papers of Fazekas, while we find 
that the reason for disagreement with the third paper- 
is their neglect of the correlation between the states of 
the nearest-neighbor bonds in counting the many-particle 
configurations. This omission of the correlation, both in 
the Gutzwiller as well as in the slave-boson mean-field 
methods, effectively turns the bond problem into the site 
problemjii^ yielding results for the ID spinless fermion 
model very similar to the ID Hubbard model. 

Gutzwiller wave Junction - The Gutzwiller variational 
wave function is written as 



(2) 



where |\l/o) is the uncorrelated wave function, g \s a, vari- 
ational parameter between and 1, and the Gutzwiller 
factor , with the bond occupancy operator B = 
^ fiifii+i, reduces the weight of the configurations where 
neighboring sites are occupied {B = 1). The uncorrelated 
wave function |^'o) can itself be taken as variational as 
welli^ The variational parameter g is determined by min- 
imizing the total energy ('H)g = (*G|?i|^'G)/(*G|*G), 
which however is difficult to evaluate analytically for 
any given |4'o) except for special cases. One then re- 
sorts to the Gutzwiller approximation for the expec- 
tation values of the one-particle operators: {0)g = 
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FIG. 1: (Color online) Initial and final configurations, and 
r^, /i — 1,...,4, that contribute to the kinetic energy term 
{c\cj). The change in the total number of bonds between the 
initial and the final configurations are indicated as well as the 
total number of configurations Kf_i{B, N, L) in each case. 



(^'g|0|^'g)/(*g|*g) ~ 7(^'o|0|*o>/(*o|*o>, where 
the "reduction factor" 7 is independent of the j'l'o) used 
in the Gutzwiller wave function. This is the Gutzwiller 
approximation. However, if one uses a |5'o) for which 7 
is exactly evaluated, then the results are fully variational 
as is the case here. 

Gutzwiller's original approach for computation of the 
reduction factor turns out to be^S exact if one uses the 
equal-weight many-body configuration for the uncorre- 
lated state: 



-1/2 



(3) 



in which case the calculation of 7 boils down to merely 
a combinatorial problem. Here iVr is the total number 
of many-body configurations. The essential complica- 
tion here is that while the site occupancies in the Hub- 
bard model are specified independently of one another, 
the bond occupancies are correlated between nearest- 
neighbor bonds, making the combinatorial problem sig- 
nificantly more involved. 

Note that, as already mentioned, since we evaluate 
the reduction factor for the equal-weight state ex- 
actly and we also use the same state in the Gutzwiller 
variational wave function Q throughout this paper, our 
results remain fully variational. We have not tested a bet- 
ter variational wave function for |^'o) in Eq. ([2]), since 
our main goal here is to obtain the Gutzwiller reduction 
factor. 

Proceeding now with the evaluation of the expectation 
values, we start with the norm 



(*g|* 



(4) 



B=Q 



where K{B, N, L), evaluated explicitly in Appendix I, is 
the number of configurations with the bond occupancy 
_B, L is the number of sites on the ID lattice (periodic 
boundary condition is assumed, so that we have a ring in- 
stead of a chain), and A^(< L) is the number of electrons. 
In the thermodynamic limit, N,L,B ^ 00 and n = N / L 
is finite, the summand is a sharply peaked function of B, 



and the sum may be replaced by the largest value de- 
termined from the condition d{g'^^K{B, N, L))/dB — 0. 
This yields the relation between the variational parame- 
ter g and the bond occupancy per site b = B / which 
reads 



6(1 -2n + b) 
[n - 6)2 ' 



or inverting this, we find 



where 



6 = n + |(.g2-l)- 



p=l- [1 + 4n(.g2 - 1)(1 - n] 



1I/2 



(5) 



(6) 



(7) 



It can be easily verified that for the uncorrelated wave 
function, g — 1^ we recover the correct result lim^^i 6 = 
. The final result for the norm is 



(n - 6)2 



In view of the fact that in the thermodynamic limit, 
I^'g) is an eigenstate of the bond occupancy operator B, 
the total energy can be evaluated as 



E{B)^ 



(^Gli^l^-c) 



(*g|*g) 



(*o|*o) 



VB, 



(9) 



where 7 = {c\cj)G/ {c\cj)o is the reduction factor of the 
kinetic energy due to correlation effects. 

Gutzwiller reduction factor 7 - To obtain the expres- 
sion for 7, we compute the expectation value 

= 1^ E + •9"'^2 + gK^ + K^), (10) 

where in the configuration jF) = |F^Ffl), F^ refers to the 
configuration of the four neighboring sites shown in Fig. 
[1] while Tfj refers to the remaining sites. Only the four 
initial and final configurations contribute to the sum and 
the Ki^s are simply the total number of configurations 
Tji in each of the four cases shown in the figure. Using 
again the thermodynamic limit to replace the sum in the 
second line of Eq. (ITUl) by its largest value and dividing 
by the norm Eq. (|5]), we find 

_ (^g|cIc,|^g) Ki+g-'K2 + gK3 + Ki 



(*g|*g) K 
(6 - bg + gn){l - 2n + b - bg + gn){n - 6) 
gn{l - n) 



(11) 
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FIG. 2: (Color online) Reduction of the bond occupancy b - 



n from the uncorrelated value (feu 



(top) and the 



inverse renormalized band mass m/m* (bottom) as a function 
of the band filling n for several values of the interaction V. 
For half filling, n — 1/2, a metal-insulator transition occurs 
for V > 8t within the Gutzwiller method, where the bond 
occupancy becomes zero and the renormalized mass becomes 
infinity. 



where the final expression has been obtained by us- 
ing the results given in the Appendix for the if 's. By 
putting g — 1 for the uncorrelated state, we readily get 
(c|cj)o — n{l — n) and after some algebra an expres- 
sion for the all-important kinetic energy reduction factor 
in the Gutzwiller theory is found. It takes a much more 
complex form than the corresponding factor for the Hub- 
bard model with the on-site Coulomb interaction, viz.. 



7 



(cJc,)g _P- 2n(l - n)(l + g){p + 1 - g) 



(clcj)o 



2n2(l-n)2(5-l)(l-fg)3 



(12) 



which is expressed in terms of the single variational pa- 
rameter g. 

The expression for the reduction factor, Eq. (fT2|) . is the 
central result of this work and it is valid for any general 
filling factor n. Three mutually conflicting expressions 
were obtained in the literature earlier^'^"^* for this reduc- 
tion factor. We find that our expression agrees with the 
second paper by Fazekas^^, which has the expression for 
7 for the half-filled case only. As explained in the sec- 
tion on the slave-boson solution below, the origin of the 
discrepancy in the literature is due to the neglect of the 
correlation between the neighboring bond occupancies, 
so that the reduction factor in Ref. 13, e. g., becomes 
identical to the case for the on-site Hubbard model. 
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FIG. 3: (Color online) Ground-state energy per electron as 
a function of electron density calculated with the correlated 
(solid line) and uncorrelated (dashed black line) wave func- 
tions. The red dots indicate the results obtained from exact 
diagonalization on a 16-site lattice with periodic boundary 
condition. The Coulomb interaction parameter used here is 
V/t = 4. 



Total energy - The total energy E{g,n) may be read- 
ily written using the Eqs. dH]), pT|) . and (O for any 
band filling n (which is the same as the electron den- 
sity). The ground-state energy is obtained from the con- 
dition dE/dg = 0. The results for the band occupancy b 
and the renormalized mass m* /m = 7uncorr/7 are shown 
in Fig. [21 where a Brinkman-Rice^^ type metal-insulator 
transition is indicated for the half- filling (n = 1/2) be- 
yond the critical value of the interaction, V > 8t. At the 
Brinkman-Rice transition point, the effective mass m* 
goes to infinity and the band occupancy goes to zero, as 
seen from Fig. [21 

The calculated total energy as a function of the band 
filling is shown in Fig. [31 where it is also compared to 
the exact results for the 16-site lattice obtained from the 
exact diagonalization using the Lanczos method. Note 
that the Gutzwiller energy forms an upper bound to the 
exact energy here, since with the equal-weight state as 
the uncorrelated state for |^'o) in Eq. the method 
is fully variational without any approximation, because 
all expectation values are evaluated exactly in this case. 
Only when we use a different uncorrelated state for |5'o), 
but at the same time still use the expression for the re- 
duction factor 7 obtained using the equal-weight state 
(Gutzwiller approximation), that the variational princi- 
ple does not hold any more. Therefore, in both the Figs. 
[3] and [51 the Gutzwiller energy forms a strict variational 
upper bound to the exact energies. 

Now returning to Fig. [31 the correlation energy, as ex- 
pected, is the largest at half filling, reducing gradually 
to zero for n = or 1. The energy for the uncorre- 
lated state was obtained by evaluating the kinetic en- 
ergy for the state j^*) = H^'^'^c^jvac), which is simply 
Eke — YTh^'^ 2tcosfc, and adding to it the interaction en- 



ergy Vb, where b 



is the bond occupancy for the 



uncorrelated state. The correlation energy, defined to be 
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FIG. 4: (Color online) Correlation energy per electron ob- 
tained from the difference between the ground-state total en- 
ergy and the same for the uncorrelated state. Correlation 
energies obtained using the ground-state energy from exact 
diagonalization (red dots) and the Gutzwiller approach (solid 
line) are both shown. 



the difference between the exact energy and the uncor- 
related energy, is plotted in Fig. ID It is clear that the 
Gutzwiller approach is able to account for most of the 
correlation energy. 

Half-filled case - The above expressions are valid for 
any general filling factor n. We now specialize to the 
half-filled case, n — 1/2, and compare with the differing 
results of the earlier authors, the resolution of which was 
in fact one of the motivations for this work. By substitut- 
ing n = 1/2 into Eqs. ([6]) and ([12]), we readily obtain the 
results for the bond occupancy and the reduction factor 



and 



& = 2-ig(l+5)- 



7 = 8g{l + 9y 



(13) 



(14) 



The total energy per site in the half-filled case, easily 
found with the help of Eqs. ([9]) and ((TT|) . is given by 



E 



Agt 



(l+.g)3 2(1 +.g) 



V. 



(15) 



Minimizing the energy with respect to the variational 
parameter dE/dg = 0, one finds that at the minimum 



(^-^+3)^/^ _i 
90^ ^175 1-^^ ' 



(16) 



if w = V/{8t) < 1, and go = 0, ii v > 1. The ground-state 
energy is plotted in Fig. [S] as a function of the strength 
of the interaction V. 

Note that for the case V — 0, go — 1/2 and the kinetic 
energy reduction factor becomes 7 = 32/27, i. e., greater 
than one, a result that has already been noted in the 
literatureJ^ At first site, this appears counterintuitive 
because one does not expect the energy of the uncorre- 
lated state to be changed by the Gutzwiller variational 
wave function when = 0, so that 7 should equal one. 
However, this need not be considered surprising since the 
equal- weight uncorrelated wave function |^'o), used to ob- 
tain the reduction factor 7, does not necessarily have the 
lowest kinetic energy, even for the case of = 0. There- 
fore the Gutzwiller variational method tries to produce a 
better ground state and ground-state energy within the 
variational degree of freedom it is given. The kinetic en- 
ergy of the correlated state could therefore in principle be 
larger than the uncorrelated trial state with the kinetic 
energy reduction factor exceeding one as is the case here. 

We illustrate this with a specific example. We consider 
the four-site, two-electron ring model. For this case, one 
can readily show for V = that the equal-weight state 
l^-o) ^ Q-^/^ ^21=1 iTfc) has the kinetic energy of -4t/3, 
while the Gutzwiller variational wave function, |^g) = 

|ri)+|r2)+<7o(|r3) + |r4)+|r5) + |r6)), where ffo = i/V2 

and the last four configurations have the bond occupancy 
B = 1, while for the first two, B = 0, has a better 
kinetic energy —\/2t. Thus the reduction factor, which 
is given by the ratio of these two energies, becomes 7 = 
(-V2t)/(-4V3) w 1.06. The exact ground-state energy 
—2t is much lower. The factor 7 is greater than one, 
although different from the infinite-lattice value of 7 = 
32/27, which it will approach as the number of sites in 
the ring is increased from four to infinity. 

Exact solution - The spinless fermion model in ID has 
been exactly solved for the half-filled case by invoking 
its equivalence with the XXZ spin Hamiltonian via the 
Jordan- Wigner transformationi^i The equivalent Hamil- 
tonian is: H = H, 



spin 



VL/4, where 



with Jx 



Sf^,S!) + jJ:^,S:, (17) 



2t, J, 



V, and L is the number of sites 



in the ID lattice. Using its exact solution obtained from 
the Bethe ansatz,~ the energy for the spinless ID model. 
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results do not change. Furthermore, a similar correlation 
is missing within the Ogawa^ treatment of the Gutzwiller 
problem, which also produces the same reduction factor 
Eq. ([T^ .^'^ This is to be compared with the equivalent re- 
duction factor from our treatment, readily obtained from 
Eqs. © and dnj: 



4 
V/2t 



8 



FIG. 5: (Color online) Exact ground-state energy per site 
for the spinless fermion Hamiltonian compared with the same 
obtained from the Gutzwiller and the slave-boson mean-field 
methods. 



Eq. ([Ij, may be written asii 



■^-2t sin 61 / (1 
J2t(i-ln2) 



tanhoj^ , 



djjj 



It sinh ( 



tanh CJTT 
— tanhn0) -|- - 



if ^ y < 
if y = 2t, 

if 2i , 
(18) 



where B = co^-^{1-^V jt) and = co&\-^(2-^V jt). 
The exact energy is plotted in Fig. [S] together with the 
Gutzwiller energy, Eqs. and (ITS)) , as well as the re- 
sults of the slave-boson mean-field theory for comparison. 

Slave-boson mean-field solution - The Gutzwiller re- 
duction factor for the half- filled case, Eq. agrees 
with the results of Fazekas^^, but not with the other two 
papers in the literature i^^'^^ The reason for disagreement 
with the last paper is due to the neglect of correlation be- 
tween the occupancies of the neighboring bonds in Ref. 
[l^ . which in effect turns the intersite interaction problem 
into an onsite problem. The slave-boson mean-field ap- 
proach neglects the same correlation as well. In order to 
understand this point, recall that in the latter method, 
auxiliary boson fields are introduced corresponding to the 
empty, singly, and doubly occupied bonds, which must 
satisfy certain constraints for the physical Hilbcrt space. 
In the mean-field theory, if these constraints are satisfied 
on each bond, without regard to the state of the neigh- 
boring bonds, it leads to the kinetic energy reduction 
factor in the slave-boson theory 



7sb 



{n - b)[Vb \/l - 2n + bf 
n{l — n) 



(19) 



which is the same result for the Hubbard model with on- 
site interaction, obtained by using either the Gutzwiller 
or the slave-boson treatment of the problem.— 12£ We have 
checked that even if we keep nearest-neighbor bond cor- 
relations within the slave-boson mean-field theory, the 



7 ^ 



{n-b)^[Vb + y/l-2n + b]^ 
n^(l — n)^ 



(20) 



where we have expressed the reduction factor in terms of 
the bond occupancy b rather than the Gutzwiller varia- 
tional parameter g. The ground-state energy per site for 
the slave-boson case is given by 



Esh = TsbCo + Vb, 



(21) 



where eo is the exact uncorrelated ground-state kinetic 
energy, eo = —2/tt at half filling, and the energy has to 
be minimized with respect to b. The results are shown in 
Fig. [S] together with the exact and the Gutzwiller results. 
2t , Wigner crystallization - Unlike the repulsive Hubbard 
model in ID at half-filling, where the ground state is 
insulatingiS for all non-zero U/t, the spinless fermion 
model does have a metal-insulator transition for the pa- 
rameter V/ {2t) — 1. The Gutzwiller approach does give a 
Brinkman-Rice type metal-insulator transition, although 
the transition point occurs at V/{8t) — 1, which is much 
higher than the exact result. The transition may be seen 
from the bond occupancy 6, which is readily calculated 
by substituting the minimum value go in the expression 
Eq. (Hg, with the result b = 6-^(2 - y^T+W/Si). The 
bond occupancy is a diminishing function of V and be- 
comes zero at and beyond the transition point, leading to 
the "Wigner crystal," where alternate sites are occupied. 

Summary - In Summary, we obtained the Gutzwiller 
solution of the spinless fermion model in ID containing 
the nearest-neighbor Coulomb interaction. Results were 
presented for a general band filling and compared with 
the exact solutions. It was shown that the Gutzwiller 
energy compares very well with the exact energy for 
intermediate coupling. Combinatorials arising in the 
Gutzwiller solution of the same problem in higher di- 
mensions appear to be quite difficult, but would be a 
worthwhile problem to study in the future. 

This research was supported by the U.S. Department of 
Energy, Office of Basic Energy Sciences, Division of Ma- 
terials Sciences and Engineering under Award No. DE- 
FG02-00ER45818. 



Appendix A: Combinatorials 

Here we outline the results for the combinatorials 
needed for the evaluation of the expectation values for 
the Gutzwiller wave function for the ID spinless fermion 
model with the nearest-neighbor interaction. The sym- 
bols Ki{B, N, L) and K{B, N, L) represent the number of 



6 



configurations for the linear chain and the ring arrange- 
ments, respectively, where N objects are arranged on a 
lattice of L sites with B bonds between them (a bond is 
defined to be present if two contiguous lattice sites are 
occupied) . 

1. Enumeration of Ki for the chain arrangement 

First of all, note the fundamental counting result that 
the number of ways to distribute n identical objects 
among k distinct containers with at least one object in 
each container is (JJZi)- Each such distribution can be 
associated with a string of n zeros and k — 1 slashes with 
each slash inserted between two zeros. The above com- 
binatorial represents the number of ways fc — 1 positions 
can be chosen for the slashes out of the n — 1 positions 
between the zeros. 

Now, we turn to the enumeration of the number of 
configurations Ki{B, N, L) for the linear chain. When 
the occupied positions form a contiguous string (with no 
intervening unoccupied positions), then B = N — I. If 
the occupied positions are separated (by unoccupied po- 
sitions) into two contiguous strings, then B = N — 2. 
In general, if C is the number of contiguous strings 
into which the occupied positions are separated, then 
B^N-C. 

For a given value of C, the number of ways to partition 
N occupied positions into C contiguous strings is (^Zi) ■ 

This is the same as (jv-s^i) = (^b ^) (since C = N-B). 

The number of ways to distribute C contiguous strings 
among L — N unoccupied positions (so that the contigu- 
ous strings are separated from each other by unoccupied 
positions) is {^~c^^) = (^n-b^)- (There are L - TV -f 1 
available locations in which we can insert the C contigu- 
ous strings of occupied positions.) Therefore, the total 
number of ways to specify N occupied positions among L 
total positions so that there are exactly B bonds is given 

by 

Note that the total number of ways to arrange the N 
objects on L sites irrespective of the number of bonds is 
simply: YJbZI Ki{B,N,L) = (^) = Nr, as it must be, 
where in the last expression iVr is the total number of 
many-body configurations. 

2. Enumeration of K for the ring arrangement 

We choose one position on the circle (it doesn't matter 
which one) and designate it as po- We now consider two 
cases. 

1. Position po is unoccupied. The number of ways 
to distribute the N occupied positions among the 



remaining positions on the circle can be obtained by 
using the formula obtained for linear arrangements 
in Eq. (|A1[) . In this case, we have a total of L — 
1 positions. Therefore, the number of admissible 
arrangements is (^"j^Is^^) = 

2. Position po is occupied. Start with a circle of N 
occupied positions, including po- We need to break 
the circle into N ^ B contiguous strings by insert- 
ing unoccupied positions among the occupied po- 
sitions. This means that we need N — B break 
points. There are N possibilities for break points. 
Therefore, there are (^^5) = (2) ways to choose 
the break points. Having chosen the N — B break 
points, there are now (^Zb~i) ways to insert the 
L — N unoccupied positions. So the total num- 
ber of circular arrangements with po occupied is 

fN\ /L-N-l\ 

The total number of circular arrangements for a speci- 
fied value of B is then given by the sum of the above two 
results. By direct computation, it can be shown that this 
sum is equal to 

W,iV,X)^^r-^V^-^Y (A2) 
^' ' ' L-N\ B J\N-Bj' ^ ' 

which is the desired result. 



3. Enumeration of the number of ring 
configurations Kij with two adjacent positions 
specified 

Here, we identify two adjacent positions on the ring 
and specify in advance whether or not each of those po- 
sitions is occupied. We designate the number of configu- 
rations by Kij, where i,j = 1 or 0, where i indicates the 
occupancy of one site and j indicates the occupancy of 
the adjacent sites. It clearly does not matter which pair 
of adjacent sites is occupied from the symmetry of the 
ring. 

1. The two positions are both occupied. The argument 
is very similar to that of case (2) in the subsec- 
tion above. Suppose positions po and pi are occu- 
pied. Now, consider a ring of N occupied positions, 
including po and pi. We need to break the ring 
into N — B contiguous strings by inserting unoccu- 
pied positions among the occupied positions. This 
means that we need N — B break "points. There are 
— 1 possibilities for the break points (since we 
can not break between po and pi). So, there are 

(jv-s) ~ (^B ^) "^^ys to choose the break points. 
Having chosen the N — B break points, there are 
now (jyZbZi) ways to insert the L — N unoccupied 
positions. Therefore, the total number of arrange- 
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merits with both po and pi occupied is 



KuiB,N,L) 
= jK{B,N,L). 



N ~1 
B - 1 



L~ N ~l 
N-B-1 



(A3) 



2. The first position is occupied and the second is not. 
Then iV — 1 occupied positions and L — N — 1 un- 
occupied positions remain to be assigned. We need 
to identify N — B — 1 break points in the occu- 
pied positions and then distribute the unoccupied 
positions. The number of ways to do this is 



Kio{B,N,L) = 



N - 1 
N~B-1 
N-l\/L-N-l 
B )\N - B-l 
N- B 



L-N-1 
N- B-l 



-KiB,N,L). 



(A4) 



3. The first position is unoccupied and the second is 
occupied. This is obviously the same as case (2), 
viz., Koi{B,N,L) = Kio{B,N,L), as dictated by 
symmetry. 

4. Neither of the two positions is occupied. This is 
similar to case (1) in subsection 2 above. The num- 
ber of arrangements is 



Koo{B,N,L) = 



N-l\fL-N-l 



B J \ N- B 
L-2N + B 



L 



-KiB,N,L). 



(A5) 



Note that when we add up these four results, we im- 
mediately see that Ku + Kio + Kqi + Kqo = K, which 
is the total number of configurations, irrespective of the 
occupations of the two chosen adjacent positions. 



4. Enumeration of the number of ring 
configurations Kij^i with four adjacent positions 
specified 

Extending the notation of the last subsection, Kijki 
denotes the number of configurations, where ijkl (each 1 



or 0) denote the occupancies of the four consecutive sites. 
Specifically, we evaluate the quantities Kiou, which we 
will need in the computation of the expectation values for 
the kinetic energy operator (c^Cfc). The arguments are 
analogous to those used in the cases above. The results 
are: 



Kioii{B,N,L) = 

B{N - B - 1 









N - 




\B- 




\N- 


B- 





-KQiiB,N,L), 



KooiiiB, N,L) = 
K,oio{B,N,L) = 
KooioiB, N,L) = 



{N - 1){L- N - ly 

N -2\/L- N -2\ 
B-l)\N-B-l)' 

B )\N - B -2)' 
N-2\/L-N-2^ 
B )[n-B-1, 



(A6) 
(A7) 
(A8) 
(A9) 



These four values are denoted as Ki , K2 , if 3 , and K4 , 
respectively, in Fig. [T] in the main body of the text. They 
can also be written as fractions of Kqi as we have shown 
explicitly for the first case and these fractions add up to 
one, so that i^ioii + -ft^ooii + i^ioio + i^ooio = Kqi, as 
it must since Kqi is the total number of configurations 
without any reference to the occupancies of the first and 
the fourth sites. 



These quantities are difficult to enumerate on the lin- 
ear chain as the circular symmetry is broken and they be- 
come dependent on the positions of the contiguous sites 
on the linear chain. Because of this, it is easier to evalu- 
ate the Gutzwiller energy on the circular ring rather than 
the linear chain. Both should of course converge to the 
same result in the thermodynamic limit. 
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